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Abstract 

The simple plant location problem is formulated as the minimization of a pseudo-Boolean 
function. This form of the problem is then transformed into a set covering problem and also into a 
weighted vertex packing problem on a graph. These formulations are compared to similar 
formulations in the literature and to the "standard” integer programming formulation. 
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INTRODUCTION 



In this paper the simple plant location problem, SPLP, is formulated as the minimization of 
a pseudo-Boolean function. This formulation is then transformed into two other discrete 
optimization problems: a set covering problem and a weighted vertex packing problem on a graph. 
These three formulations of the SPLP are compared to similar formulations that have appeared in 
the literature and the differences are discussed. 

The SPLP is described as follows. Let P = {i: i = 1, . . , p} be the index set of potential 
locations for plants (or plant sites) in some space such as a network or the plane. If a plant is 
opened at location i, a fixed cost fj is incurred. Let D = { j: j = 1, . . , d) be the index set of 
customers. The unit cost of transportation between customer j and plant location i is given by cjj 
for each i 6 P and je D. Each customer has a demand which must be met by the opened plants 
and is assumed to be one unit. If a customer's demand is different from one unit, it is scaled to one 
and the transport costs are scaled accordingly. Assume further that the capacity for each opened 
plant is sufficient to meet the demand of all customers (hence the alternative name "uncapacitated 
plant location problem"). The SPLP is then to choose a subset of locations from P at which plants 
are opened and to specify the transportation between opened plants and all customers so as to meet 
customer demand and to minimize the total fixed cost plus the total transportation cost. 

As an example, consider the network in Figure 1 where there is a customer at each node and 
a plant site at each of the nodes 1, 3 and 4. Each edge number is the transportation cost accross that 
edge. For each plant site i and customer j, the transportation cost Cjj is the minimum cost over all 
paths between i and j. Table 1 gives the fixed cost fj and the transportation costs c- for each i e P 
and j e D. 
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Figure 1: A tree network with customer and plant sites at nodes 
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Table 1: Costs data for the SPLP of Figure 1 



The approaches discussed here are not limited to a SPLP defined on a tree network, nor even to a 
network. The approaches apply to any SPLP defined by a matrix of nonnegative cy values and 
nonnegative fixed costs fj. 

In the literature there are several recent surveys of location problems that include a 
discussion of the SPLP. The extensive surveys by Comuejols, Nemhauser and Wosley [3] and by 
Krarup and Pruzan [8] are devoted almost entirely to the SPLP. Each includes a thorough 
discussion of the problem's origins and several of its formulations. References [2,3] consider the 
"uncapacitated facility location problem" which is equivalent to the SPLP, but maximizes the total 
revenue from satisfying customer demands minus the cost of opening plants at the chosen sites. 
The formulations discussed below may be applied directly to the uncapacitated facility location 
problem with only minor modifications. 



A BOOLEAN FORMULATION 

We formulate the SPLP directly as the minimzation of a pseudo-Boolean function. This 
formulation follows naturally from the well known property that in some optimal solution, each 
customer j will receive its entire unit of demand from one open plant, namely a plant with minimum 
transportation cost to customer j. 

For each i e P, define the variable yj to be 1 if a plant is opened at location i, and 0 
otherwise. Let yj = 1 - yj be the complement of y- v For each customer j, let j(.) be a 
permutation of the location indices i e P so that the transportation costs from each j(i) to j are in 
nondecreasing order: 

c j(l)j - c j(2)j - - c j(p)j' ^ 

Observe that in meeting the demand of customer j, the cost c j(i)j is incurred iff a plant is 
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open at j(l), i.e., yj^) = 1; or the cost Cj( 2 )j is incurred iff a plant is open at j(2) and no plant is 
open at j(l), i.e., yj( l)Yj(2) = or g enera F the cost c j(k)j i ncurre d iff a plant is open 
at j(k) and no plant is open at locations with smaller cost: j(l),j(2), .. ,j(k-l), i.e., yj( l)yj(2) * • • 
yjCk-l)yj(k) = L Finally, if no plant is open, i.e., yj(i)yj(2) • • • y j(p) =1 ’ the demand of 
customer j is not met. This corresponds to an infeasible solution and is avoided by including in the 
problem a penalty term Mjyj ( i ) yj (2) • • • yj(p) w f' ere Mj is a large cost. Suitable values for the 
Mj are discussed subsequently. 

To simplify the notation, let Flyj^^ denote the product yj(l)yj(2) * • • yj(k) f° reac h 
j € D and k = 1, . . , p. Furthermore, let rTyj(O) = 1 in subsequent expressions. Then the 
transportation cost to customer j from the set of open plants is given by 

T j vS^iAo-i/jd) + M i n *j (P ) • <2) 

In the example problem, for customer 2, we have 



T 2 = ly 3 + 2y 3yi + 4y 1 y 3 y 4 + M 2 y 1 y 3 y 4 . 

The total fixed cost incurred by the open plants is given by the sum of the fjyj. Then, a 
pseudo- Boolean function F (see reference [7]) of the variables yj, y 2 , . . ,yp is defined as the total 
fixed cost plus the sum of the Tj's for each j in D: 

F( y l’ • • • ’ y p > ” 2r f i y i + C j (i )j ^ y j (i- 1 ) y j (i) + M j ny j(p) 1 



The SPLP may be written as 

PI: min { F(yj, . . ,y p ) : yj = 0,l, ieP}. 

A formulation equivalent to PI is reported by Hammer [6] and also included in the survey [8]. 

In a specific instance of the SPLP, the function F can generally be simplified by adding 
together those terms with identical variables. Two customers j and j' will have terms with identical 
variables, say Flyj(k-l)yj(k) Tj> and FIyj'(k-l)yj'(k) Tj'’ whenever the sequence of plant 
sites, ordered by j(i) and j'(i), are identical for i = k, k+1, . . ,p for some k < p. This ability to 
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aggregate like terms is a desirable feature of the pseudo-Boolean function F that is not found in 
most other formulations of the SPLP. In particular, since each penalty term ffy^pj is identical for 
each customer j, all these may be added together to yield one penalty term. Let M denote the sum 
of the Mj over j e D. Then the function F for the example problem, simplified by adding together 
like terms, is given as follows: 

F(y],y 3 ,y 4 ) = 5y 1 + 7y 3 + 5y 4 + y^ + 3y 3 y j + 5y 4 y 3 + 1 ly 3 y iy 4 + 7y 4 y 3 y! +My 1 y 3 y 4 . 



PENALTY VALUES 



It is convenient to rewrite the function F, using the common penalty term M, as: 

P dp 

F( y r • • ■ . V - IVi + I |l c j(i)j n Wj<i) + Mny j(p) • 



The following discussion shows how the value of F depends on the penalty term M. Observe that 
a set of variable values (y j , . . . , yp) is feasible to PI if and only if Flyj(p) = 0. Furthermore, if 
(y j , . • . , yp) is an infeasible solution, then F(yj, . . . , yp) = M, the penalty value. Consider a 
minimum feasible solution (y*j, . . . , y*p), and let z* = F(y*i, . . . , y*p) denote the optimal 
objective function value. 

The first claim is that with penalty value M < z*, each minimum solution to F is infeasible. 
This follows since each infeasible solution has the objective function value M < z* < F(y j, 
yp) for each feasible solution (yj, . . . , yp). Alternatively, with a penalty cost of M > z* each 
minimum solution to F has optimal value z*. This follows since for each feasible solution 



(yj, . . . , yp), z* < F(y], . . . , yp), and for each infeasible solution (y’j, . . . ,y'p), z* < M = 
F(y’j, . . . , y'p). This same argument shows that with penalty cost M > z*, each minimum 
solution to F is feasible. 

Thus for all M > z*, the optimal objective function value of PI is z*, and any upper bound 
on z* will suffice for a value of M. Values for the penalty coefficients Mj may be obtained from 
any feasible solution (y* j , . . . , y*p). Define 



f* 



P 

= 

itl 1 



and c*j = min ( c j(i)j ny* j(i _ 1)y * 



jO) 



i = L 



). 



that is, f* is the total fixed cost of this solution and c*j is the minimum transport cost to customer j 
for this solution. Then each penalty cost Mj may be chosen as any value greater than c*j + P/d. 
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A SET COVERING FORMULATION 



Problem PI is next transformed into a set covering problem. This formulation has fewer 
variables and constraints than a set covering formulation of the SPLP by Kolen [9], while it retains 
the combinatorial properties of his formulation. It is also related to a covering formulation for the 
uncapacitated facility location problem reported in reference [2] but obtained by a different 
approach. 

For each variable yj^ which appears in uncomplemented form in expression (2) substitute 
yj (i ) = 1 - yjQ and simplify. Expresion (2) then becomes, for customer j, 

p-1 

T j ~ C jO)j + i?/ C K i+1 )j ’ C j(i)P n ^j‘(i) + ^ M j * C j(P)P n ^j(P)' ^ 

In the example problem, for customer 2, we have 



T' 2 = ! + y 3 + 2y 1 y 3 + (M 2 - 4)y 1 y 3 y 4 . 



Then the function F has an equivalent form, denoted by F’, which is the sum of the Fixed costs plus 
the sum of the T j for each j in D: 



P d 




The function F has a natural interpretation. For each customer j e D, the cost c j(i )j is a 
constant that is incurred by any solution to the SPLP. That is, the minimum transport cost to 
customer j, c j( l)j, must always be incurred. This cost may be zero if a plant site and a customer 
location coincide. Each term (Cj(i + j)j - c j(i)j) is the incremental transportation cost incurred by 
customer j if plants at j( 1), j(2), . . ,j(i) are not opened. In this case, the transportation cost to j 
is at least Cj(i+i)j, and if j(i+l) is open the cost is exactly Cj(j+i)j. 

As above, the function F’ may be simplified by adding together like terms. For the 
example problem this yields: 
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F (yi>Y3’y4) = 2 + 5yj + 6y 3 + 4y 4 + y\ + 2y 3 + 4y 4 + 7 y^ + 2y 3 y 4 + (M - 18)y !y 3 y 4 . 



A further reduction is possible for F' between pairs of complementary variables yj and yj 
using the identity yj + y; = 1. Any two terms of the form ayj+ byj may be rewritten as follows: 

ayj+ by j = { (a-b)yj + b, if a > b; or Oyj + (b-a)yj + a, if a < b } . (4) 

In the second case, an uncomplemented variable yj with zero cost may be given the value one, and 
all terms containing the complementay variable yj will have value zero and may be eliminated from 
F. For the example problem the reduction (4) yields: 

F '(yi>y 3 >y4) = 9 + 4y\ + 4y 3 + 0y 4 + 7y 1 y 3 + 2y 3 y 4 + (M - 18)y 1 y 3 y 4 . 

Here y 4 is retained in the function F' in order to illustrate the next transformation. 

Next, the SPLP is transformed by replacing each term Fly^ in F by a variable znj(i)’ 
where Flj(i) = (j(l)j(2). . .j(i)}, and by adjoining the constraints 

z rij(i) ^ n yj(i) i e P, j e D. (5) 

But each constraint of (4) is equivalent, for zero-one variables and their complements, to the 
constraint 

i 

/ y VI , + z >1 

nj«) 



These transformations, applied to the general expression for F', yield the set covering 
problem: 



p 


d p-1 


d 


P2: min j/.y. 
i=l 


+ j?J ^ c K i+1 )j" c W 


1 z nJ(i) + ^ M j ’ c j(p)? z nj(p)^ + j?j c j( 1 )j 


s. t. 


i 

^ y j(k) + Z nj(i) ” 1 


i E P, j E D 




yi ’ z nj(i) = ^>i 


ie P, je D 
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Problem P2 will have at most d(p-l) + p + 1 variables and d(p-l) + 1 constraints. 

When F is simplified by adding together like terms and by applying the reduction (4), the 
number of variables znj(I) num ber of constraints in P2 will be decreased accordingly. For 

the example problem, with F' simplified and M' = (M - 18), problem P2 becomes: 



min 9 + 4yj + 4 y 3 + 0y^ + 7 zj 3 + 2 Z 34 + M'zj 34 




s. t. yj + y 3 


+ z 13 


> 1 


y3 + y4 


+ z 34 


> 1 


yi + y3 + y4 


+ z 134 


> 1 



all variables zero - one. 

The above example suggests further simplifications that may apply to specific problems. 

Any yj variable with a zero cost coefficient may be set to one in any optimal solution. The 
corresponding column and all rows with a one in that column may be eliminated. In the example 
problem, setting y 4 = 1 and eliminating columns and rows leaves one constraint: yj + y 3 + zj 3 > 1 , 
from which it is easily seen that the two alternative optimal solutions are to open plants at sites 
{1,3} or {1,4} with a cost of 13. Thus the set covering formulation may be interpreted as 
preprocessing of the problem data that may specify some plant sites to be opened a priori . 

The well known set covering 'reduction rules' [10] will generally apply to problem P2. For 
example, since the column corresponding to a variable yj always dominates a column 
corresponding to a variable Zfij'(k) f° r 1 € ^j(k), the column of zj-jj^ can be eliminated if the cost 
coefficient of yj is not greater than that of zj-jj^). In particular, since the coefficient of yj is less 
than M, the variable zpy^p) and its column can always be eliminated. 

Row reductions may follow. If the variable zj-jj^ has been eliminated from a row by a 
column reduction, this row may now be dominated by some other row. In this case, the dominating 
row may be eliminated. 

Problem P2 is quite similar to Kolen's set covering formulation of the SPLP [9] with the 
following differences. First, Kolen's formulation does not allow the simplifications we get from 
adding together like terms. Each such simplification eliminates at least one variable and constraint. 
His formulation has a penalty variable zp-jj^p-j and its corresponding constraint for every j in D, 
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whereas these terms are combined into one term and constraint in problem P2. These observations 
imply the constraint matrix of problem P2 is a proper subset of the constraint matrix of Kolen's 
formulation. The example problem is taken from Kolen [9]. His set covering formulation of the 
example problem requires 17 constraints and 20 variables. 

The set covering reduction rules also have more limited application to Kolen's formulation. 
In general, these rules would eliminate the columns corresponding to the variables z nj(p) f° r eac h j 
in D, but since other terms are not agregated, fewer columns may be eliminated. For the example 
problem, only the rows and columns corresponding to the z nj(p) variables could be eliminated. 
Furthermore, in Kolen's formulation, no yj variables have zero coefficients so that no plants can be 
set open as in the example problem. 

Kolen showed that the constraint matrix of his set covering formulation was totally balanced 
[9] when the transportation costs were weighted distances in a tree network and the nodes served 
both as customer and plant locations. Since the constraint matrix of P2 is a submatrix of the matrix 
of Kolen's formulation, and since a submatrix of a totally balanced matrix is totally balanced, this 
same property holds for the formulation P2. 

Problem P2 is closely related to a formulation reported in reference [2] for the uncapacitated 
facility location problem and obtained by a canonical reduction of the matrix of cy values. This 
formulation also aggregates terms where possible and so reduces the number of variables and 
constraints. 



A WEIGHTED DOMINATION PROBLEM 

The set covering problem P2 is equivalent to a weighted domination problem on a bipartite 
graph. To see this equivalence, problem P2 is expanded by adding redundant variables and 
constraints so that the expanded matrix of constraint coefficients is the vertex-vertex incidence 
matrix of a bipartite graph. 

Consider the function F', which has been simplified by adding together like terms and 
reduced by expression (4). Let S be the subset of indices i in P corresponding to those terms in F' 
that contain only y;. For each i in S, yj represents the aggregation of one or more variables yj ( q ) 
over those indices j such that j(l) = i. To obtain problem P2 each variable yj was replaced by the 
variable zj and the constraint yj + Zj > 1 is included. 

A new problem, called P2', is obtained by expanding problem P2 as follows. For each i in 
P-S, add the variable Zj to problem P2 with a zero cost in the objective function and add the 
constraint yj + zj > 1. Problem P2' now has a variable zj and a constraint yj + z\ > 1 for each i in 
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P. This version of P2' is equivalent to P2 since any feasible solution to P2 may be amended by Zj 
= 1 for i in P - S to obtain a feasible solution to P2' with no change in the objective function. 
Conversely, deleting zj = 1 from any feasible solution to P2' gives a feasible solution to P2 with 
no change in the objective function. 

Let m be the number of variables z nj(i) * n including the Zj added above. Note that m is 
also the number of constraints in P2'. The coefficient matrix of P2' may be written as (I m ,A), 
where I m is an mxm identity matrix with a column corresponding to each variable z nj(i)’ an d A is 
an mxp matrix with a column corresponding to each variable yj. Observe that A contains the pxp 
identity matrix Ip as a submatrix since for each i in P, there is a constraint zj + yj >1. 

Next, the set of constraints with coeficient matrix (A 1 , Ip), where A 1 is the transpose of A, 
is added to P2'. Then the coefficient matrix of P2' becomes: 

c _ ( M 

v At 'p ; 

For each i in P, the i^ constraint of (A 1 , Ip) contains Zj and yj (since A contains Ip) and is therefore 
satisfied by any solution that satisfies the constraint zj + yj > 1 in (I m ,A). Thus the constraints of 
(A 1 , I p ) are redundant and problem P2' is equivalent to P2. 

Since the matrix C is symmetric with ones on the diagonal, it is a vertex-vertex adjacency 
matrix of some graph G. The graph G is seen to be bipartite with one set of vertices corresponding 
to variables yj and the other to the variables zj-jj^. Each vertex zjqj(i ) ls a djacent to those vertices 
yj where i is in the set I~Ij(i). Each vertex is weighted by the cost coefficient of the respective 
variable. Therefore, the set covering problem P2' is a weighted domination problem on the 
bipartite graph G. Figure 2 shows the bipartite graph of the example problem. 
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Figure 2: Bipartite Graph of the Example Problem 



A VERTEX PACKING FORMULATION 

The third formulation obtained for the SPLP is that of a weighted vertex packing problem 
on a graph. For each customer j and its permuation j(i) for i = 1 , . . . , p, we have the following 
Boolean identity: 



J, nii j(i-i) y j(i) + n h( P ) = 1 



Solving for ITyj^, substituting into expression (2) and simplifying yields 



T " 

J 



P 

J/ C j(i)j ’ M j )ny j(i-l) y j(i) 



+ M. 
J 



(7) 



for each j in D. For the example problem and customer 2, 



T"2 = (1-M2)y3 + (2-M2)y3yi + (4-M2)yiy3y4 + M2. 



Next substitute Vj= 1 - yj in each term fjyj of F. Making these substitutions in the function F 
gives an equivalent form: 



F(y r • • 





dp P 

- y y (M. - c.,. v )ny v - ,,y v -,. + m + Vf- . 
pi^i J JWJ J j(i-irj(!) pi 1 



Define the function G(y j, . . . , y p ) = - F(yj, . . . , y p ), which is a pseudo-Boolean function in 
posiform. Then problem PI may be written in the equivalent form as: 



P3: max {G(yj, . . . , y p ): yj = 0,1, ieP}. 

The function G for the example problem, with Mj = 5 and without simplifications, is: 

G(yj, Y3,y4) = 5y\ + 6y 3 +4y 4 + 5yj + 4y jy 3 + yjy 3y 4 + 4y 3 + 3y 3 y ] + ¥1^4 + 
5 y3 + 4y3 yi + 2y 1V3V4 + 5y 4 + 2y 4 y 3 + j^j^yi + 4y4 + 3V4V3 + 2V3V4V j - 40. 
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Hammer and Rudeanu [7], and Ebenegger, Hammer, and de Werra [5] have shown that the 
problem of maximizing a pseudo— Boolean function given in a posiform may be transformed into a 
maximum weighted vertex packing problem on a graph. In reference [5], the graph corresponding 
to the posiform pseudo- Boolean function is called a "conflict graph". The conflict graph has a 
vertex for each term in the function, and an edge connecting two vertices if the corresponding terms 
contain at least one complementary pair of variables. That is, if yj is in some term of G and its 
complement yj is in some other term, there is an edge between the two vertices. The vertex weights 
of the graph are the coefficients of the terms in G. We call the conflict graph of the function G 
derived from a SPLP a "plant location graph" (PLG). The plant location graph for the expample 
problem, without any simplifications is shown in figure 3. 

It is shown in reference [5] that the maximum value of G is equal to the value of a maximum 
weighted vertex packing of the conflict graph. We denote the problem of determining a maximum 
weighted vertex packing in the plant location graph of the function G as problem P4. 

The SPLP has been formulated as a weighted vertex packing problem by Cho, Johnson 
Padberg, and Rao [1] and by Cornuejols and Thizzy [4], Our formulation is closely related theirs 
but with some interesting differences. An obvious similarity is the highly structured form of the 
plant location graph. For each j in D, the set of vertices corresponding to the terms in T'j form a 
clique (’rows' in Figure 3). For each i, a star is determined with the vertex corresponding to yj at 
the center and an adjacent vertex corresponding to the term containing yj in each expression T"j 
('columns' in Figure 3). 

Several differences are observed. First, the plant location graph of P4 contains all the edges 
of the graph constructed in references [1,4] but has the additional edges (between cliques) 
corresponding to complementary variables. For example, there is an edge between y^ and yjj^ 
in the example problem. These additional edges strengthen the vertex packing formulation. 

Second, if no simplifications are made in G by adding together like terms, the plant location 
graph has the same number (dp+p) of vertices as does the graph constructed in [1,4]. However, 
each simplification of G by adding like terms reduces the number of vertices in the plant location 
graph. This reduction in vertices is not possible for the graph constructed in references [1,4]. 

Third, the reduction given by expression (4) may be applied to the function G. For each i 
consider the terms (Mj - c j(i)j)yj(i) + fjyh suc h that j(i) = i. For the function G, the penalty cost 
Mj may be chosen sufficiently large so that Mj - cj^j > fj. Then these two terms reduce to (Mj - 
c j(i)j ' fj)yj(i) + fj which eliminates the vertiex corresponding to to yj from the plant location graph. 
With these reductions, the plant location graph has at most (dp) vertices. The conflict graph 
resulting from the function G, after all simplifications and reductions, is called the reduced plant 
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location graph. 

For the example problem without simplifications in the function G, the plant location graph 
contains 18 vertices and 73 edges. However applying all possible simplifications and reductions, 
the function G becomes: 

G(yi, Y 3 » y 4 > = 3y 3 + 5y 4 + 4yjy 3 + 7y 1 y 3 + 5y 3 y 4 + 4y L y 3 y 4 + 3y!y 3 y 4 - 25. 

Figure 4 shows the reduced plant location graph for the function G as simplified above. 

Finally, for all SPLP with p plant sites and d customers, the vertex and edge sets of the 
graphs constructed in references [1,4] are identical. Only the vertex weights change. However, a 
plant location graph may differ from one SPLP to another with the same sets P and D since the 
vertex set and the edge set depend directly on the costs cy and fj. In particular, the plant location 
graph for a SPLP arising from a tree network may differ from a SPLP arising from a general 
network of the same size. This raises the following question corresponding to Kolen's result for 
the SPLP on a tree network: Does the plant location graph arising from a SPLP on a tree network 
have some special struture? 




RELATIONSHIP TO STANDARD INTEGER FORMULATIONS OF THE SPLP 

Problems PI and P3 are equivalent to the well known integer programming formulation of 
the SPLP, often called the "strong formulation" [3] which is given below. The variables yj are the 
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same as defined above. Define Xjj to be a binary variable indicating whether customer j is served 
by a plant at i (xy = 1) or not (xjj = 0). Then the strong formulation may be written as follows: 



P5: minSX, c ij x ij +Y f i>i ® 

1=1 J=1 J J 1=1 

P 

s. t. X x ii = 1 j 6 D (9) 

i=l J 

xjj < yj j e D, i € P (10) 

xjj = 0,1 yj = 0,1 (11) 



Given next is a relationship between the variables xjj and the Boolean expressions in the yj 
variables that were defined earlier. For any set of values for the yj variables, xjj = 1 if and only if a 
plant is open at site i and the transportation cost Cjj between customer j and site i is minimum over 
all other open plant sites. That is, for each i e P and j e D with the j(.) defined by expression (1) 
and j(k) = i, 



x ij — yj(i)yj(2) * * * yj(k-i)yj(k) - n yj(k-i)jyj(k). ( 12 ) 



Expression (12) is substituted for each Xjj in problem P5. Thus each constraint of (10) 
becomes 

n Vj ( k-i)jyj(k) S Vi 



which is always satisfied for yj equal to either 0 or 1. With these substitutions, the constraints (10) 
can be discarded since they are always satisfied. 

Substituting expression (12) for each xjj into constraint (9) yields 



d 

k ?, nSi j(k-i) y j(k) " ’ j£D (13) 



Comparing (13) to the Boolean identity (6) implies that each constraint (13) is equivalent to 
requiring that klyj(p) = 0 for each j in D. Thus the constraints (9) may be eliminated from P5 if 
each term flyj(p) is forced to be zero in an optimal solution. This is accomplished by adding to the 
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objective function the terms MjETyj^pj where Mj is a large penalty cost, for each j in D. 

Finally, the cy may be reindexed by the permutations j(.) to yield cjyy, and expression 
(12) is substituted for each xy in the objective function. With these changes, problem P5 is written 
as 



dp p 

min £ w Wjo) + Mn W + 

s. t. yj = 0,1 

which is precisely problem PI . 

Since problems PI and P3 are equivalent, it follows that P5 is equivalent to P3. However, 
it is interesting to note that P3 may be interpreted as a Lagrangian relaxation of P5. Expression 
(12) is again substituted for xy in the objective function (8) and in each constraint (9) and (10) of 
problem P4. The constraints (10) are deleted as above since they are always satisfied. For the 
resulting problem, consider a Lagrangian relaxation with respect to the constraints (9) and with 
multipliers Mj. This yields the following problem. 

p d p d p 

min ilijli c j(i>j ni j(i)j y jao + + £, M i (l ' L ny j(k-i) y j(k>) 

Yi = 0,1. 

Problem P3 is obtained by substituting yj = 1 - yj in each term fjyp collecting terms, and 
changing to a maximization problem by multiplying through by -1. 



15 




Figure 3: Conflict graph of G without simplifications 
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